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0\ : Abstract 

Tunneling between two 2D electron gases in a weak magnetic field is of 

in 

resonance character, and involves a long lifetime excitonic state of an electron 
and hole uniformly spread over cyclotron orbits. The tunneling gap is linear 

oo 

ly^ ' in the field, in agreement with the experiment, and is anomalously sensitive to 

o ■ 

' the electron density mismatch in the wells. The spatial coherence of tunneling 

in 



o 



along the orbit can be probed by magnetic field parallel to the plane, which 



' produces an Aharonov-Bohm phase of the tunneling amplitude, and leads to 

a" 

an oscillatory field dependence of the current. 



Electrons in GaAs quantum wells have very high mobility, and at low temperature form an 
almost ideal Fermi liquid. Strong Coulomb interaction and sensitivity to external magnetic 



^ ' fields make the physics of this system rich and interesting. Recently, several phenomena 



have been discovered in tunneling experiments in quantum wells, including: a tunneling 
gap induced by a magnetic field resonance peaks of conductivity near zero bias 0; and 
excitonic effects |^ . Typically, one has two wells containing 2D Fermi liquids with an electron 
density of the order of lO^^cm"^, separated by an oxide barrier of few tens nanometers thick. 
It is characteristic that the barrier thickness is very uniform, so that tunneling is coherent in 
lateral dimensions. This results in conservation of both momentum and energy. Basically, 
at a zero magnetic field, only tunneling between identical plane wave states can occur, which 
restricts the phase space of final states and leads to an unusual I — V curve with a sharp 
peak near zero bias p|. Finite width of the peak is determined by elastic scattering. One 



can probe spatial coherence by applying a magnetic field parallel to the barrier. Having 
practically no effect on the dynamics in the plane, the vector potential of the parallel field 
shifts electron momenta in one plane, which creates a mismatch of the Fermi surfaces of 
different planes. The effect on the tunneling rate can be easily accounted for by the free 
electron picture p. 

A magnetic field applied perpendicular to the plane makes the situation more interesting, 
especially in high fields, when the system is in the Quantum Hall state [Q . Compared to the 
I—V curve at zero field, the current peak is shifted away from zero bias to some finite voltage. 
Also, the peak broadens as the field increases. The current is almost entirely suppressed 
below the lower edge of the peak, called a "tunneling gap." Recently the tunneling gap has 
been intensively studied because it is believed that it can be used to probe the QH state. 
The gap depends on the field hnearly at weak field (z/ ^ 1) 0], and saturates at higher field 
(z/ ^ 1) [||]. To summarize, one can say that the I — V curve can be interpreted in terms of 
a resonant tunneling mechanism, involving some intermediate state with the lifetime given 
by inverse peak width, and of the work needed to create this state corresponding to the gap. 

The gap at high field is quite well understood its energy scale is of the order of 

e^/ea, where a = n~^/^ is the interparticle distance, and e is the dielectric constant. Un- 
til recently, the low field gap has received less attention. The only available theory is by 
Aleiner, Baranger, and Glazman |^, who developed a hydrodynamical picture by treating 
the system as an ideal compressible conducting liquid. They wrote down classical elec- 
trodynamics equations in terms of charge and current densities, and derived the tunneling 
gap A = {huc/ 1') In^ue^ /ehvp). At constant electron density, they predict quadratic depen- 
dence of the gap on the magnetic field, which is different from the linear dependence found 
experimentally 

The reason for the disagreement, in our opinion, is that in a clean metal, such that uJcT S> 
1, one can use classical electrodynamics only on a scale much bigger than the cyclotron 
radius Rc = vp/uJc- However, we will argue that the important scale of the problem is of 
the order of Rc, and thus one has to have a Fermi liquid in magnetic field. On this scale, 
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the state formed at tunneling has a non-trivial spatial structure, which makes the physics 
very different from that of the gap in high fields. Also, we will find that at weak field a large 
simplification occurs, because the problem is semiclassical, and one can use the classical 
Fermi liquid equation to describe the dynamics in terms of Fermi surface fluctuations. 

Another point is that in two dimensions the energy and momentum conservation pre- 
scribes that all quantum numbers of the final and initial tunneling states coincide. In a 
magnetic field, this implies that the radius of an electron orbit as well as the guide center 
of the orbit are conserved at tunneling. In a weak field, this results in spatial coherence of 
tunneling over a large distance of the order of Rc, and leads to interesting effects. 
Summary of results. Our goal is to explain the linear field dependence of the gap in a 
weak magnetic field, and to propose experiments that will reveal spatial structure of the 
intermediate tunneling state. Although at the end we are going to do a rigorous Fermi 
liquid calculation, it is instructive to begin with a semiclassical picture of electron states 
localized near classical cyclotron orbits and weakly interacting with each other. At tunneling, 
an electron hops from the orbit of radius Rc in one layer to the identical orbit in the other 
layer, and leaves a hole on the first orbit. This creates an electrostatic configuration of two 
oppositely charged rings of radius Rc separated by the barrier of width d -C Rc- The energy 
of this charge distribution is 

enRc Lb 



where e is the dielectric constant, and the magnetic length Ib = yhc/eB characterizes the 
ring's "thickness." Basically, we are saying that, in order to transfer an electron, one has to 
charge the "two ring capacitor," and its charging energy e^/2C constitutes the tunneling gap. 
The result (|1]) holds for < d, i.e., fields which are not too low, and fall in the experimental 
range [Q. The gap A dependence on magnetic field is nearly linear, since the log term is 
roughly constant. The dependence on the barrier width is in agreement with the excitonic 
picture By the order of magnitude the gap (|ip agrees with the experiment which 

raises the question of why there is no gap suppression due to Coulomb screening. Basically, 



the reason is that he charge istribution is locahzed in a very thin ring, of the thickness Ib 
comparable to the screening length, which makes the screeing effects not too dramatic: they 
simply change the log in Eq.(|l]) by a constant of the order of one. 

That the tunneling is coherent along the cyclotron orbit can be easily verified by applying 
a magnetic field parallel to the barrier, in addition to the perpendicular field. The parallel 
field flux "captured" between electron and hole trajectories will give an Aharonov-Bohm 
phase to the tunneling amplitude, and make it proportional to 

J J e^RcdB^^icos9i-cose2)^Q^^Q^ _ (2) 

Thus the current will oscillate as a square of the Bessel function: 

m = Ji /(v)b„=o . (3) 

For an ideal two dimensional system, the current dependence on B\\ factors out because 
a parallel field does not affect the motion in plane. For real wells of finite width, the 
factorization (j^) should still be a good approximation at small 5^, when the cyclotron 
radius is big compared to the well width. However, since the parallel field will squeeze the 
states in the wells, and effectively increase the barrier width, the tunneling rate may aquire 
an additional non-oscillatory suppression factor. 

Also, it is clear from what has been said, that the gap will be very sensitive to any 
asymmetry between the wells. For example, if there is a small mismatch 5n of densities in 
the wells, the radii of the electron and hole orbits will become different, 5Rc — {Rc/2n)6n, 
which will reduce mutual capacitance of the orbits, and thus raise the electrostatic energy. 
The gap will be enhanced by the order of its magnitude at 6Rc — Ib, which corresponds to 
the density mismatch 6n/n ^ Ib/ Rc anomalously small at weak field. 
Fermi liquid calculation. To study the screening effects, we will use the Landau equation 

(ai-D(l + F))MP,r,t) = , (4) 

where D = Vp ■ Vr + uocmvp x Vp contains the Lorentz force term. In Eq. we have omitted 
a collision integral, which is legitimate in the semiclassical limit hujc <^ Ep. We will find a 

4 



stationary solution that fully accounts for the screening of the electron and hole charges. The 
resulting distribution consists of two parts: singular and smooth. The singular distribution 
is localized near the electron and hole cyclotron orbits, within a short distance of the order 
of magnetic length. The smooth part is spread over a distance larger than the cyclotron 
radius. It is worthwhile to make a comparison with the standard picture of a particle in 
a Fermi liquid, which consists of a quasiparticle excitation combined with a background 
density fluctuation, and to draw a relation with the singular and smooth densities. 

For simplicity, let us assume that F corresponds to a pure density-density interaction: 

(FHa(P, r) = E U^P^^ - ^')^^/3(p', r') , (5) 

where a, f3 = R, L label wells, and Uapir) is the Coulomb potential. It will be straightforward 
to modify the calculation for a more general Fermi liquid interaction. One can write 5n(p, r) 
corresponding to the electron and hole being uniformly spread over cyclotron orbits: 








where r(t), p(t) is the classical circular trajectory with |p| = pp, and with the period 
T = 2ti/uJc. Semi classically, the density (|^) corresponds to a state of the lowest Landau 
level available for tunneling. At tunneling, initially, the states (^) are created, and then, 
over the time scale ~ "-u"^, they relax to a stationary state, whose energy gives the tunneling 
gap. 

In order to find the system response to the appearance of the electron and hole, we have 
to solve the equation: 

{dt - D{1 + F)) 5n = 5{t) 5n(°) . (7) 

Formally, the solution of Eq. (|^) can be written as a "ladder" : 

Sn= (^G + GDFG + GDFGDFG + ...) , (8) 

where G = {dt — D)^^ is readily evaluated in the Fourier representation: 
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in{0p—9pi)+i\iX{p—p')/mLJc 

G(k,p,p')=E TT-. ^ • (9) 

The density satisfies the equation with F = 0. The first term of the ladder 
gives 5n^^\p,r) 9{t), the singular part of the density. The rest is the smooth "background" 
density. 

In terms of 6n the gap is given by 

A = J2 (P, r)F"^5n^(p', r', t ^ oo) . (10) 

p,r,a,p' ,r' ,f3 

By summing up the ladder (H) in a standard Fermi liquid fashion, one gets 

A = E Sn%M^^iSn^^^,ipX^o , (H) 

k,p,p' 

(l + F,)J-^,k = Afl+^E^^^-^-,kl , (12) 

Here J„ are Bessel functions. At u; ^ it turns into = Fk (l + Fk{l — JlikRc))^ . By 
plugging it in Eq.(|ll]) together with the Fourier transform of density, / c^^P^^^k r,(l)(p) = 
± Jo{kRc), we finally get 

J J',{kR,)e'U, kdk 

where Uk = 27r(l — e~^'^)/ek is the 2D Fourier transform of the difference of the in-plane and 
interplane Coulomb interaction, and v = mj'nV? is the compressibility. 
Discussion of Eg. (fT^j. We can evaluate the contributions of two regions in the fc-space, 
kRc ^ 1 and kRc <C 1, by replacing Jo{x) with its asymptotic expressions: 

d 

(i) AkR^^i = — - — hucln—; (14) 

nehvp Tg 

AkR^«i = ^ In - . 15 
mvp Tg 

According to the experiment 0], we now assume that Ib < < d <^ Rc, where Ib is 
magnetic length, and is the Coulomb screening length. The contribution ([M]) coincides 
with the estimate (|1]), except for the log term suppressed by screening. When the result 



(0) is compared to the experiment , one has to be careful because the screening length is 
probably set by the well width, rather than by the electrons' compressibility. Because this 
may alter the log term, we can claim only an agreement with the experimentally measured 
A by the order of magnitude. 

It is interesting to note that the contribution ( [T^ ) coincides with the gap found by 
hydrodynamics theory the difference in the log term being due the effect of screening 
by another plane. The two terms (|T^) , (p!5l) account for contributuions of the two parts 
of the charge density, the singular and the smooth. Not surprisingly, the singular density 
dominates in the energy, which one can understand by comparing charging energy of a 
two-ring capacitor with that of a parallel plate capacitor, both of the size Re- 
Bosonization calculation: a sketch. Now we proceed with an accurate derivation of tunnel- 
ing current. In the weak field, huo^ ^ Ep, the electron motion is semiclassical. Therefore, 
instead of doing full many-body theory, we can write down the action that corresponds to 
the classical Landau Fermi liquid equation, and use it to study dynamics in imaginary time. 
By that, we will determine optimal path in imaginary time (instanton) whose action will 
give the exponent of the tunneling rate. 

The advantage of this approach, besides making a clear relation with classical theory, is 
the possibility of including effects of a parallel magnetic field and of a scattering by disorder 
in a natural fashion. 

We use a bosonized Fermi liquid picture to write electron operators in terms of 

Bose fields: 

^«'''W = 7!^/'"""*"'"'°'""'°- '''' 

Here the unit vector n = p/p^ labels points of the 2D Fermi-surface, and tq is a cutoff of 
the order of the inverse bandwidth that appears in the bosonization formalism. The integral 
/ ...(in means / ...c?6'/27r, where 9 is the polar angle in the p— plane. 
The bosonized imaginary time action for 0(n, r, t) is 

S = ^J(f)D(idt-D{l + F))(j)dndhdt. (17) 



Here (j) = {(pji, (pL), and F is the integral operator (H). 

Note that u = -D0(n, r, t) has a clear meaning of the displacement of the Fermi-surface at 
the point n in the normal direction, i.e., along n. The saddle points 0(n, r,t) of this action 
satisfy Eq.(^). According to Haldane in the context of the Fermi liquid bosonization, 
the operator D has an interpretation of a covariant derivative. 

To get the tunneling current, one has to find the amplitude 

K{v, r) = (T.^+(r, zr)^^(0, 0)M^, zr)^+(0, 0)) . (18) 

Then, one continues K to the real time axes from the upper and the lower halfplane. (Fol- 
lowing we denote these functions as K''^^\v,t). By the standard formalism |jl2|, one 



can express the tunneling current as 

/ = e|topIm J cPx dt ($r,te~'"^* + <l>_r, V^*) , (19) 

where "I'r,* = K^{Y,t) — K^{r,t). To find K, one has to evaluate the functional integral of 
the product of the exponentials ( |16|) with the weight . One notes that the integral is 
gaussian, and writes: 

K{r,t) = ^ I e-^"o."icinocini . (20) 

Here the action 

^no,n, = {JD-^ (idt - D{1 + F)) J) , (21) 

where (...) = / ...(in(i^R(it, and 

J = 5(r)5(t)(5(n-no) - -5(r-R)5(t-r)5(n-ni). (22) 

The source term J describes an electron injected at t = 0, r = at the point no of the Fermi 
surface, and then removed at t = r, r = R from a different point rii. 

Calculation of the action ( PT| ) involves inverting the operator id — D{\ \ F\ which has 
been done above (see (||), (|])), and then applying it to J. The latter procedure requires 
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a lot of attention, because of the singular nature of the operators in (pT]). A complete 
discussion of this operator, including its relation to classical trajectories, and regularization 



technique will be presented elsewhere The result of this calculation, however, is very 
simple: I{V) ~ 6{eV — A)sign{V) , where A is given by Eg. ([T3|) . In other words, the rigorous 
calculation confirms the semiclassical result for the gap, and predicts an infinitely sharp 
resonance peak at eV = A. 

To include scattering by disorder, we modify the action (|T7|) by adding a "collision 
integral" term: 

TTl P doJ 

S=— 0-..(n, r)Diy0^(n, r) signer— dn dh. (23) 
An J In 

Here W is the collision integral operator 

Wu (n) = - - -(M(n) - / u{n) dn) , (24) 

Ts J 

taken in the so called r— approximation. It is straightforward to verify that the new action 
leads to the classical Fermi liquid kinetic equation with the scattering term. Let us note that 
in the imaginary time representation the action ( P5D is non-local, because it is dissipative. 
The Gaussian functional integral over again gives (^Up, where now 



5, 



_ At^ 



{JD-^ {tdt - D{1 + F)-W) V) . (25) 



In the absence of a magnetic field, by evaluating the action (^), then plugging it in 
Eqs.(|20),(O), and doing the integral [O], we get 



' s 

The peak of conductivity dl/dV is located at zero bias, and has finite width due to elastic 
scattering, and Lorentzian tails, in agreement with the experiment MM- 



A perpendicular magnetic field creates the tunneling gap. Calculation shows that 
the I — V curve is related to the zero-field curve (B^) in a simple way: 
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Io{V - A/e) at V > A ; 



at -A<V < A] 



(27) 



loiy + Aje) at V <-A 



where A is given by Eq.(|I|). So, at finite the peaks of loiV) are simply shifted by A/e 
away from zero bias. 

This formula correctly describes the width of the gap and the peak shape, however it 
doesn't hold deep inside the gap, near y = 0, because our simple approach to scattering 
breaks down at t ^ r^. In this limit the charge spreading regime changes to diffusive, and 
one has to deal with a classical electrodynamis problem |jl3|. 

To include a parallel field, we modify the amplitude K{y, t) by introducing in pO|) the 
Aharonov-Bohm phase given by the parallel field fiux: 



where y is the component of the radius vector r perpendicular to the field B\\ . The expression 



tunneling current, which makes the current an oscillating function of B\\. 
Conclusion. We interprete the I — V curve in terms of a resonance tunneling mechanism 
that involves an excitonic state, whose structure in the weak field limit can be studied by 
the Fermi liquid theory. We calculate the tunneling gap, and find linear dependence on 
the magnetic field, in agreement with the experiment. The structure of the excitonic state 
manifests itself in the gap high sensitivity to any asymmetry between the wells, like density 
mismatch. Also, it can be probed by a magnetic field parallel to the barrier, by observing 
oscillatory dependence of the current on the field. 

We study the problem by means of the Fermi liquid bosonization theory, and obtain a 
complete I — V curve that includes the effects of the perpendicular and parallel fields, as 
well as of elastic scattering. 



K{v,t) 



(28) 



(|I^) then produces 




This gives a prefactor to the 
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